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Abstract 

We study the hadroproduction of a CP-even or CP-odd neutral Higgs boson in 
association with a Z boson in the minimal supersymmetric extension of the standard 
model (MSSM) We include the contributions from quark-antiquark annihilation at 
the tree level and those from gluon-gluon fusion, which proceeds via quark and 
squark loops, and list compact analytic results. We quantitatively analyze the 
hadronic cross sections at the CERN Large Hadron Collider assuming a favorable 
supergravity-inspired MSSM scenario. 
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1 Introduction 



The search for Higgs bosons is among the prime tasks of the CERN Large Hadron Collider 
(LHC) [1]. While the standard model (SM) of elementary-particle physics contains one 
complex Higgs doublet, from which one neutral CP-even Higgs boson H emerges in the 
physical particle spectrum after the spontaneous breakdown of the electroweak symmetry, 
the Higgs sector of the minimal supersymmetric extension of the SM (MSSM) consists of a 
two-Higgs-doublet model (2HDM) and accommodates a quintet of physical Higgs bosons: 
the neutral CP-even h° and H° bosons, the neutral CP-odd A boson, and the charged 
P^-boson pair. At the tree level, the MSSM Higgs sector has two free parameters, which 
are usually taken to be the mass m^o of the A boson and the ratio tan /3 = v 2 /v i of the 
vacuum expectation values of the two Higgs doublets. 

In the following, we focus our attention on the h° ,H°, and A bosons, which we 
collectively denote by 0. A recent discussion of P ± -boson production at the LHC may 
be found in Refs. [2-4] and the references cited therein. The dominant source of bosons 
is their single production by gg fusion, gg — > 0, which is mediated by heavy-quark [5] 
and squark [6] loops. Another, less important mechanism of single 0-boson production is 
bb — > [7]. The bosons thus produced have essentially zero transverse momentum (pr)- 
In order for the bosons to obtain finite px, they need to be produced in association with 
one or more other particles or jets (J). In leading order (LO), j0 associated production 
proceeds through the partonic subprocesses gg — > g<p, gq — > g0, and qq — > g<p, which 
again involve heavy-quark [8,9] and squark [9,10] loops. Alternatively, the bosons can 
be produced, with interesting rates, in association with (i) a dijet via intermediate-boson 
fusion, qq' — > qq'V*V* — > qq'(p, where q and q' stand for any light flavor of quark or 
antiquark, V = W ± , Z, and virtual particles are marked by an asterisk [11]; (ii) a quark- 
antiquark pair of heavy flavor Q = t,b via gg, qq — > QQ<f) [12]; (iii) an intermediate boson 
via qq' -> W ± <p [13-16], qq Z<p [13,14,16-21], and gg Z(j) [17-19,22-25]; or (iv) 
another, possibly different boson via qq — > 0102 [16,26,27] and gg — > 0102 [27-29]. 
Note that, due to the absence of A°VV couplings at the tree level, qq' — > qq'V*V* — > 
qq'cj) and the Drell-Yan processes qq' — > W ±if — > W ± (f) and qq — >■ Z* — >■ Z<p are not 
possible for = A . The partonic suprocesses gg — > Zcf) and gg — > 0102 are mediated by 
heavy-quark [17-19,22-25,27-29] and squark loops [18,27,29]. Comprehensive reviews of 
quantum corrections to Higgs-boson production within the SM and MSSM may be found 
in Refs. [30,31], respectively. 

In this paper, we revisit Z<p associated hadropro duct ion via gluon fusion in the MSSM. 
In the SM case, mutual agreement between three independent calculations [17,22,25] has 
been established, and compact formulae for the partonic cross section are available [17]. 
In the MSSM, the status is much less advanced and, perhaps, somewhat unsatisfactory. 
As for gg — > Zip with = h°,H°, there exists only one analysis so far [19], which has 
not yet been verified by other authors. The analytic expressions presented in Ref. [19] 
are rather complicated; they involve 14 form factors. In order to translate the SM results 
[17,22,25] to the case of gg — > Z(j> with = h°, H° in the MSSM, it is not sufficient to 
adjust the HZZ and Hqq couplings; it is also necessary to include certain quark triangle 
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diagrams with an A boson in the s channel [see Fig. 1(a)]. On the other hand, the 
squark loop contributions to these two MSSM processes vanish [19] for reasons explained 
below. As for gg — > ZA°, the quark loop contributions were first studied on the basis 
of a numerical evaluation [23], which was recently employed for phenomenological signal- 
versus-background analyses taking into account the subsequent Z — > l + l~ and A — > bb 
decays [20,24]. An independent analysis, including also the squark loops, was reported 
in Ref. [18], which does not contain an analytic expression for the partonic cross section 
either. Unfortunately, comparisons between Refs. [18,20,23,24] are not discussed in these 
papers. Recently, the helicity amplitudes of gg — > Z<ft with = h°, H°, A were analyzed 
in Ref. [32] with regard to the asymptotic helicity conservation property of supersymmetry, 
and their real and imaginary parts were graphically presented at a specific scattering angle 
as functions of the center-of-mass (cm.) energy. 

With the Higgs hunt at the LHC being in full swing, it is an urgent matter to con- 
solidate our knowledge of Z<f> associated hadroproduction via gluon fusion in the MSSM, 
which is the motivation of this paper. Specifically, we present compact analytic expres- 
sions, also for the partonic cross sections of qq, bb — > Z<p [18-21], and perform a detailed 
numerical analysis using up-to-date input. 

The importance of 66-initiated subprocesses for Higgs-boson production has been var- 
iously emphasized in the literature, in particular, in connection with the final states [7], 
0102 [27], H + H- [33], and W ± H^ [2,34]. These subprocesses receive contributions from 
Feynman diagrams involving 6-quark Yukawa couplings, which are generally strong for 
large values of tan /3. (The btH~ and ibH + couplings are also strong for small values of 
tan ft.) If the two final-state particles couple to a Z boson (or photon), as is the case for the 
final states h°A°, H°A°, Zh°, ZH°, and H + H~, then there are additional contributions 
from Drell-Yan-type diagrams, which are already present for the light flavors q = u,d, s, c. 
However, diagrams of the latter type are absent for the final states h°h°, h°H°, H°H°, 
A A , ZA°, and W ± H T , which can still be produced through bb annihilation. 

As for bb annihilation, it should be noted that the treatment of bottom as an active 
flavor inside the colliding hadrons leads to an effective description, which comprises contri- 
butions from the higher-order subprocesses gb — > Z<pb, gb — >■ Z<f>b, and gg — > Z<pbb. If all 
these subprocesses are to be explicitly included along with bb —¥ Z<fi, then it is necessary 
to employ a judiciously subtracted 6-quark PDF in order to avoid double counting [7,35]. 
The evaluation of bb — > Zip with an unsubtracted 6-quark PDF is expected to slightly 
overestimate the true cross section [7,35]. For simplicity, we shall nevertheless adopt this 
effective approach in our analysis, keeping in mind that a QCD-correction factor below 
unity is to be applied. In fact, such a behavior has recently been observed for bb — > ZA° 
[21]. 

In order to reduce the number of unknown supersymmetric input parameters, we 
adopt a scenario where the MSSM is embedded in a grand unified theory (GUT) involving 
supergravity (SUGRA) [36]. The MSSM thus constrained is characterized by the following 
parameters at the GUT scale, which come in addition to tan/3 and m A o: the universal 
scalar mass tuq, the universal gaugino mass the trilinear Higgs-sfermion coupling A, 

the bilinear Higgs coupling B, and the Higgs-higgsino mass parameter /i. Notice that m^o 
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is then not an independent parameter anymore, but it is fixed through the renormalization 
group equation. The number of parameters can be further reduced by making additional 
assumptions. Unification of the r-lepton and 6-quark Yukawa couplings at the GUT scale 
leads to a correlation between m t and tan /3. Furthermore, if the electroweak symmetry 
is broken radiatively, then B and /x are determined up to the sign of \i. Finally, it turns 
out that the MSSM parameters are nearly independent of the value of A, as long as 
\A\ ^500 GeV at the GUT scale. 

This paper is organized as follows. In Sec. 2, we list the LO cross sections of qq — > 
Z(f), including the Yukawa-enhanced contributions for q = b, and those of gg — > Z(f>, 
including both quark and squark loop contributions, in the MSSM. The relevant quark and 
squark loop form factors are relegated to Appendix B. In Sec. 3, we present quantitative 
predictions for the inclusive cross sections of pp — > Z<p + X at the LHC adopting a 
favorable SUGRA-inspired MSSM scenario. Sec. 4 contains our conclusions. For the 
reader's convenience, the relevant Feynman rules are summarized in Appendix A. 

2 Analytic Results 

In this section, we present the LO cross sections of the partonic subprocesses qq — > Zip 
and gg ->■ Z<p, where = h°, H°, A , in the MSSM. We work in the parton model of QCD 
with rif — 5 active quark flavors q — u,d, s, c, b, which we take to be massless. However, 
we retain the 6-quark Yukawa couplings at their finite values, in order not to suppress 
possibly sizable contributions. We adopt the MSSM Feynman rules from Ref. [37]. The 
couplings of the Z and bosons to quarks, vz qq , az qq , and g^q, are given in Eq. (5) of 
Ref. [33] and Eq. (A3) of Ref. [27], respectively. As for the (f>ZZ couplings, ghPzz and 
9h°zz are given by Eq. (A. 3) in Appendix A, while the A°ZZ coupling vanishes at tree 
level. The h°A°Z and H°A°Z couplings, gh»A°z and gH°A°z, m ay be found in Eq. (A2) 
of Ref. [27]. For each quark flavor q there is a corresponding squark flavor q, which 
comes in two mass eigenstates i = 1,2. The masses m qi of the squarks and their trilinear 
couplings to the bosons, g^q^, are listed in Eqs. (A. 5), (A. 7), and (A. 8) and in Table 1 
of Ref. [38] \ Eq. (A. 2) of Ref! [33], and Eq. (A4) of Ref. [27], respectively. 

Considering the generic partonic subprocess ab — > Z<p, we denote the four-momenta 
of the incoming partons, a and b, and the outgoing Z and bosons by p a , pb, pz, and p^, 
respectively, and define the partonic Mandelstam variables as s = (p a +Pb) 2 , t = (p a —Pz) 2 , 
and u — (pb — Pz) 2 - The on-shell conditions read p 2 a = p\ = 0, p\ = m 2 z = z, and 
p\ = mi = h. Four-momentum conservation implies that s + t + u = z + h. Furthermore, 
we have sp\ = tu — zh = N, where px is the absolute value of transverse momentum 
common to the Z and (f) bosons in the cm. frame. 

The tree- level diagrams for bb — > Z(j> with = h°,H° and = A are depicted 
in Fig. 1(a) and (b), respectively. As already mentioned above, there is no Drell-Yan 
diagram in Fig. 1(b) because of the absence of a A°ZZ coupling at the tree level. The 
differential cross sections for the first class of partonic subprocesses may be generically 

: In Ref. [38], m qi and g<j>q iqj are called Mg a and Vg ab /g 7 respectively. 
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written as 

do- 
lt 



{bb -+ Z<j>) = ^ [(2z+p 2 T ) g\ zz (v 2 Zbb + a 2 Zbb ) \V z (s)\ 2 + \\P\ 2 

g<t>bbazbb Re P + g\ bb (v 2 Zbb T + + a 2 Zbb T_ 



- 4spr ( 7 + 
\t u 



where Gp is Fermi's constant, c w = mw/mz is the cosine of the weak mixing angle, 
A = s 2 + z 2 + h 2 - 2(sz + zh + hs), and 



P = 9<i>A°z9AObbPAo(s), 



Here, 



T± = 2 ± 2 + 2p 2 T 
Px(s) 



t u) tu 



s — m 2 x + im x Tx 



(2) 



(3) 



is the propagator function of particle X, with mass mx and total decay width Fx- For 
the second class of partonic subprocesses, we have 



dcr /,- „ , n \ G\c^ UJ z 



^ (bb -+ ZA° 
dt V 



A|,S'|- - ( j + i ) ,/ vu^izuJXvS 



3tts 

(4 bb T + + a 2 Zbb T_) 



(4) 



where 



5* = gh°A°z9h°bbPh ( s ) + 9H0A0z9H f >bbPH0 (s) ■ (5) 

As for Zh° and ZH° production, there are also sizable contributions from qq annihilation 
via a virtual Z boson for the quarks of the first and second generations, q — u,d, s, c, 
whose Yukawa couplings are negligibly small. The corresponding Drell-Yan cross sections 
are obtained from Eq. (1) by putting P = T± = and substituting b — > q. The resulting 
expression agrees with Eq. (2.8) of Ref. [17], appropriate for qq — > ZH in the SM, after 
adjusting the HZZ coupling. The full tree- level cross sections are then obtained by 
complementing the 66-initiated cross sections of Eq. (1) with the Drell-Yan cross sections 
for q — u,d, s, c. 

The non- vanishing one-loop diagrams pertinent to gg — > Z<p, with = h°, H° and 
6 = A are depicted in Figs. 2(a) and (b), respectively. As already mentioned in the 
Introduction, the presence of the quark triangle diagrams involving an s-channel A°-boson 
exchange in Fig. 2(a) represents a qualitatively new feature of the MSSM as compared to 
the SM. Furthermore, similarly to Fig. 1(b), quark triangle diagrams with an s-channel 
Z-boson exchange do not appear in Fig. 2(b). In the following, we refer to a squark loop 
diagram involving an s-channel propagator as a triangle diagram. The residual squark 
loop diagrams are regarded to be of box type. The squark triangle and box diagrams for 
gg — >■ Z(j) with <fi = h°, H° vanish, and so do the squark box diagrams for gg — > ZA°. 
This may be understood as follows, (i) The g gMj , gggq^, g g zq iqj , and g ZM . couplings 
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are symmetric in i and j, while the SU ^ coupling is antisymmetric [39]. Thus, squark 
loops connecting gluons and Z bosons with an odd number of A bosons vanish upon 
summation over % and j. (ii) The g g q iqj and gzq^ couplings are linear in the squark four- 
momenta, while the gggq^, ggZqiqj, and g^qj couplings are momentum independent [39]. 
Thus, a squark loop connecting gluons, Z bosons, and <fi bosons vanishes upon adding its 
counterpart with the loop-momentum flows reversed if the total number of gluons and Z 
bosons is odd. 

As in Refs. [3,4], we express the quark and squark loop contributions in terms of 
helicity amplitudes. We label the helicity states of the two gluons and the Z boson in 
the partonic cm. frame by A a = ±1/2, X b = ±1/2, and X z = 0,±1. We first consider 
gg — > Z(f) with = h°,H°. The helicity amplitudes of the quark triangle contribution 
read 



M 



A 



A a A o 



= -2W-(A a ±A 6 )£ 



aZqqg^ZzVzis) (F A (s, 77lj) + 2) 



—gAOqqg ( />AOzVAo(s)F A (s,mj) 

771/7 



(6) 



The quark triangle form factor, F A , is given in Eq. (B.7). As for the quark box contribu- 
tion, all twelve helicity combinations contribute. Due to Bose symmetry, they are related 
by 



■M Xa \ b x z 



(-l) Xz Ml XaXz (u,t), 
(t,u) = M\_ Xh _ x (t,u). 



Keeping X z — ±1 generic, we thus only need to specify four expressions. These read 



(7) 



zX 
8 



9HQ a Zqq m q [ F l+ + 
9Hq a Zqqm q K~ ~ (t Ii) 



<2N 



M+ +Xz = -4W — J2 g<j>qqazq q m q Fl + - (t u) 



+-A 2 



'2N 



4\/ 9<t>qi a Zqq m q F \_ + (t -B- U, X Z ~X Z ) 



(8) 



The quark box form factors, F X X J^ , are listed in Eq. (B.8). For the reasons explained 



above, we have M A 



Aa A(,Az 



■M-\ a x b x z = for the squark- induced helicity amplitudes. 



We now turn to gg — > ZA°. The helicity amplitudes of the quark and squark triangle 
contributions read 



M 



A a A b 



-8J-(1 + X a X b )J2m q (g h °A°Z9h°qq'Ph°(s) + gH°A°ZgH°qqV H °{s)) F A (s, mfj , 
V z a 
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2\/-(l + A a A 6 ) (gh°AOZ9h°qm'P h o(s) + gH°AOzgH°q i qi'PH°(s)) F A (s, 771?) , 



(9) 



(10) 



respectively. The quark and squark triangle form factors, Fa and F&, may be found in 
Eq. (B.9). Again, the helicity amplitudes of the quark box contribution satisfy the Bose 
symmetry relations of Eq. (7). We find 



M n ++0 



9A° qq az qq m q + (t u) 

i 

J2 9A° qq a Zqq m q + (t u) 



I2N 



- UI — J2 9A° qq azq q m q - (t u) 



I2N 



- i \l — 9A0 qq a Zqq m q [F+_ - (t u, X z ->■ -X z ) 



(11) 



The quark box form factors, F^jj , are presented in Eq. (B.10). We recall that Mx a \ b \ z 
0. 

The differential cross section of gg — > Z(p is then given by 



256(4tt) 3 s 2 



(12) 



where a s (fj, r ) is the strong-coupling constant at renormalization scale /i r . Due to Bose 
symmetry, the right-hand side of Eq. (12) is symmetric in t and u. 

The differential cross section of gg — > ZH in the SM is obtained from Eqs. (6)-(8) 
and (12), with <fi = h°, by replacing h° — > H, adjusting the h°ZZ and h°qq couplings, 
and discarding the contribution due to A°-boson exchange. In this way, we recover the 
result of Ref. [17], which is expressed in terms of Lorentz-invariant form factors rather 
than helicity amplitudes. 

The kinematics of the inclusive reaction AB — > Z<f) + X, where A and B are colliding 
hadrons, is described in Sec. II of Ref. [2]. Its double-differential cross section cPa/dydpT, 
where y and pt are the rapidity and transverse momentum of the Z boson in the cm. 
system of the hadronic collision, may be evaluated from Eq. (2.1) of Ref. [2]. 



3 Phenomenological Implications 

We are now in a position to explore the phenomenological implications of our results. 
The SM input parameters for our numerical analysis are taken to be Gj? = 1.16637 x 
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1(T 5 GeV- 2 , m w = 80.399 GeV, m z = 91.1876 GeV, m t = 172.0 GeV , and m b (m b ) = 
4.19 GeV [40]. We adopt the LO proton PDF set CTEQ6L1 [41]. We evaluate a s (fi r ) 
and mb(fir) from the LO formulas, which may be found, e.g., in Eqs. (23) and (24) 
of Ref. [42], respectively, with rif — 5 quark flavors and asymptotic scale parameter 
Aqcd = 165 MeV [41]. We identify the renormalization and factorization scales with the 
Z<p invariant mass y/s. We vary tan/3 and m^o in the ranges 3 < tan/3 < 32 m t /mi> 
and 180 GeV < m^o < 1 TeV, respectively. As for the GUT parameters, we choose 
mi/2 = 150 GeV, A = 0, and \x < 0, and tune m so as to be consistent with the 
desired value of m^o. All other MSSM parameters are then determined according to the 
SUGRA-inspired scenario as implemented in the program package SUSPECT [43]. We 
do not impose the unification of the r-lepton and 6-quark Yukawa couplings at the GUT 
scale, which would just constrain the allowed tan /3 range without any visible effect on 
the results for these values of tan/3. We exclude solutions which do not comply with 
the present experimental lower mass bounds of the sfermions, charginos, neutralinos, and 
Higgs bosons [40]. 

We now study the fully integrated cross sections of pp — > Z<p + X at the LHC, with 
cm. energy = 14 TeV. Fi gures 3-5 refer to the cases <fi = h°, H°, A , respectively. 
In part (a) of each figure, the dependence is studied for tan/3 = 3 and 30 while, in 
part (b), the tan/3 dependence is studied for m A o = 300 and 600 GeV. We note that the 
SUGRA-inspired MSSM with our choice of input parameters does not permit tan /3 and 
m A o to be simultaneously small, due to the experimental lower bound on m h o [40]. This 
explains why the curves for tan/3 = 3 in Figs. 3-5 (a) only start at m A o « 280 GeV, while 
those for tan/3 = 30 already start at m A o « 180 GeV. 

In Figs. 3 and 4, which refer to </> = h°, H°, respectively, the total gg-annihilation con- 
tributions (dashed lines), corresponding to the coherent superposition of Drell-Yan and 
Yukawa-enhanced amplitudes, and the g g-fusion contributions (solid lines), which arise 
only from quark loops, are presented separately. For a comparison with future experimen- 
tal data, they should be added. For comparison, also the pure Drell-Yan contributions 
(dotted lines) are shown. As for </> = h°, we observe from Fig. 3 that the contribution due 
to qq annihilation is almost exhausted by the Drell-Yan process and greatly exceeds the 
one due to gg fusion, by a factor of 3-5. The qq- annihilation contribution falls off by a 
factor of two as m h o runs from 82 GeV to 115 GeV and feebly depends on tan /3, except for 
the appreciable rise towards the lower edge of the considered tan (3 range. The gg-fusion 
contribution feebly depends on m h o, m A o, and tan/3. The situation is very different for 
(j) = H°, as is obvious from Fig. 4. Here, bb annihilation is generally far more important 
than the Drell-Yan process, except for m^o = 300 GeV and tan/3 = 3, where the latter 
gets close. The contribution due to bb annihilation monotonically increases with tan (3, 
while the one due to the Drell-Yan process decreases. Furthermore, gg fusion competes 
with qq annihilation and even dominates for tan f3 <,7. 

As for = A , the 66-annihilation contribution (dashed lines) and the total gg-fusion 
contribution (solid lines), corresponding to the coherent superposition of quark and squark 
loop amplitudes, are presented separately in Fig. 5. For comparison, also the gg-fusion 
contribution due to quark loops only (dotted lines) is shown. As in the case of = H°, 
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gg fusion competes with bb annihilation and even dominates for tan [5 <^ 7. Again, the bb- 
annihilation contribution monotonically increases with tan/3. The bulk of the gg-fusion 
contribution is due to the quark loops, especially at low values of m^o. 

Finally, we compare our results with the literature. As already mentioned in Sec. 2, we 
recover the well-known SM result [17], for = H, by taking the SM limit of our results for 
<fi = h° in Eqs. (6) and (8). The contribution due to A°-boson exchange in Eq. (6), which 
is not probed in the SM limit, agrees with the analogous contribution to gg — > W~H + 
given in Eq. (1) of Ref. [3] after appropriately adjusting the masses and couplings. On 
the other hand, the residual terms in the latter equation, which arise from the exchanges 
of h° and H° bosons, coincide with Eq. (9) after substituting the appropriate masses and 
couplings. Similarly, by adjusting masses and couplings in Eq. (10), we reproduce Eq. (2.3) 
in Ref. [4], which gives the squark triangle contribution to gg — > W~H + . In Ref. [23], 
numerical results for the cross section of pp — > ZA° via quark-loop-mediated gg fusion 
were presented. Adopting the input parameters and proton PDF set specified in that 
reference, we nicely reproduce the separate contributions due triangle and box diagrams 
shown in Fig. 4 therein, while we fail to agree with their superposition. Furthermore, 
we find reasonable agreement with the cross section of pp — > ZA° + X via gg fusion 
represented graphically for different scenarios in Figs. 6 and 7 of Ref. [21] adopting the 
respective inputs from there. 

4 Conclusions 

We analytically calculated the cross sections of the partonic subprocesses qq — > Z<p and 
gg Z(j), where = h°, H°, A , to LO in the MSSM. We included the Drell-Yan and 
Yukawa-enhanced contributions to qq annihilation (see Fig. 1) and the quark and squark 
loop contributions to gg fusion (see Fig. 2). We presented these results as helicity ampli- 
tudes expressed in terms of standard scalar one-loop integrals. 

We then quantitatively investigated the inclusive cross sections of pp — > Z<fi + X at the 
LHC with y/S = 14 GeV adopting a favorable SUGRA-inspired MSSM scenario, varying 
the input parameters m^o and tan (3. Our results are presented in Figs. 3-5. The total 
cross section for <fi = h° is typically of order 1 pb, while those for <fi = H°, A are of 
order 100 fb (10 fb) for m^o = 300 GeV (600 GeV). Assuming design luminosity, L = 
10 34 cm~ 2 s _1 , a cross section of 1 pb corresponds to 10 5 events per year and experiment 
at the LHC (see Table I of Ref. [44]). 
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A Feynman rules 

In this appendix, we collect the Feynman rules used in this paper. The Feynman rules 
for the Zqq vertices are ig^(vz qq — CLzqqls), with g = e/s w , e being the proton charge, 

s l = 1 - 4> and 

Iq — 2s w Q q I q . . 

V Z qq= , a Z qq = - — , (A.l) 

where I g = ±1/2 and Q q = 2/3,-1/3 are the weak hypercharge and electric charge of 
quark q, respectively. The Feynman rules for the 4>qq ((f) — h°, H°) and A°qq vertices are 
igg^qq and gg A o qq lb, respectively, with 

m t cos a m t sin a m t cot j3 

9h°tt = — 7, '■ — qi 9H°tt = -7, : — qi 9a°u = 



2m w sin j3 "" 2m w sin j3 " ' 2m w 

nib sin a m b cos a m b tan (3 

9h%b = 7, ni 9H%b = -7, o> 9A%b = 7, > ( A - 2 ) 

2mw cos p 2mw cos p 2mw 

where a is the mixing angle that rotates the weak CP-even Higgs eigenstates into the 
mass eigenstates h° and H°. The Feynman rules for the cf>ZZ vertices are igg<f>zzg >iu , with 

9h°zz = -— sin(a - ff), gmzz = — cos(a - f}). (A.3) 
c w c w 

The Feynman rules for the <f)A°Z vertices are gg^A°z{p + p'Y ■, where p is the incoming 
four-momentum of the <j) boson, p 1 is the outgoing four-momentum of the A boson, and 

cos(o; — j3) sin(a — j3) , 
gh°Aoz = — - , gmAoz = — - • (A.4) 

The Feynman rules for the (pqiqj vertices are igg^qj, with 

/ 9h°hii 9h ( >t 1 t 2 | 
\9h0hh 9hH 2 t 2 J 

( m z sin(q+/3) (l?-s w Qt) _ m\ cos a m t (p, sma+A t cos a) \ _ T 

Cw mivsin/3 2m w s\nP J (_A/fM 

mt sin a+At cos a) mz sin(a+P)s w Qt tn t cos o I \ / ' 

2mw sin c w mw sin / 

/ fl , /i°&i&i fl7i°&ib 2 ] 
\9h°hbi 9}fib 2 b 2 J 

( m z sin(a+/3)(j3-s^,Q t ,) mg sin a r n b (p cos a+A b sin a) 



' mw cos ft 2mw cos P ^ ] (_A/f^ 



m b (fi cos a+A b sin a) m z sm(a+0)s w Q b . m\ sin a 

\ 2mw cos Cw mw cos P 



r 
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f 9H°hti 9H°hi 2 \ 

\9h°Mi 9H°hh J 

(m z cos(a+/3)(/ t 3 -sJ,Qt) ^ m^sina m t {ii cos a- A t sin a) 

c w m w sm/3 2m w sin | (_/Vl*) 

m t (p. cos ot- At sing) m z cos(g+/3)s^,Qt _ m| sina I V / ' 

2mvy sin /? c™ sin /3 . 

9H% 2 bi 9H°b 2 b 2 J 

m z cos(o!+/3)(jg-sg J Q t ,j m| cos a mi,(^ sin a- A b cos a) 




m w cos/3 2m w cosf3 ( J^/[ h \ (A 5) 

m z cos(Q+/3)s?„Q ft mucosa V / ' v ' 7 



mi, (/i sin a— A;, cos a) cos(a+/3)5g J Q b _ m ft cos a 

2mw cos /3 c m ni\y cos /3 



A4* = COS l Sm ?" (A.6) 
V — sm t/g cos / 



are the squark mixing matrices, with 9q being the squark mixing angles. 



B Quark and squark loop form factors 

In this appendix, we express the quark and squark triangle and box form factors, Fa, 
and i 7 !^, for <fi = h°,H° and = A , in terms of the standard scalar three- and four-point 
functions, which we abbreviate as Q* fe (c) = C (a,b,c,m 2 ,m 2 ,ml) and D^f(e,f) = 

D (a, b, c, d, e, /, m?, m 2 , m|, , respectively. The definitions of the latter may be found 
in Eq. (5) of Ref . [3]. 

The quark triangle form factor for = h°,H° reads 



F A (s,m 2 q )=4m 2 q C™ q (s). 
The quark box form factors for <j> = h°,H° read 



1 ++ 



2s(t + u)C™ q (s) + 2[t + u+- 



TV I t + u + - I +2m^A 



(* " W^it) + (t- h)C™(t)\ 
D%£(t, u) - 4 (szh + m 2 q \) D h q ^(s, t), 



(h — z — s) 
N 



(t - u) \s{t + u)CZ(s) - \C h *{s) - 2m 2 ND^(t, u) 



+ 2(t + u)(t- z) 
2(t - h) 



1 + 



qqq\ 

t(t-u)(h 



mi 
z — s 



+ 



N 

— (h — z — s 



N(t + u) 

z(u 2 -t 2 -X) + (t + u)(t 2 - zh)] C™(t) 
t 2 -zh 



2st- 



N 



+ 4m (t - u) 



jjhzOO, 
U qqqq 



s,t), 



(B.7) 
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F 1 



(t-u) 



z — h — s 



Vx 



"A; 



Dqqqq ( S ? t) 



2(h-u) 



+ - — ' I X z \f\~ + t-u + 



2N 



(h - t)CZ(t) + (z - u)Cf(u) 



K- = £ + v^A - Xz(t - «)) C™(s) - | (VA + A z (f - «)) C7j» 

21 7 ^ + 3t) -2N + (u-t)(t-z) + X z (t - s - z)Vx) C™ q (t) 



VXN 

+ 2 ^_T^ (3u(s - z) + */i - 2s(/i - 2m) - A z (/i - u)>/A) C; 



zO 
999 



+ 



[t (A + 8z/i - Us - 2(t + u){z + h) 



y/XN 

+ X z (-2h + 3t + u - 2z)Vl) - l6m 2 q N 



jjhzOO, 

+ \ (-Vx - ^ + x z (t - «)) ^-°(t, «). 



M) 



The quark and squark triangle form factors for = A read 



F A ( S ,mj)=2+(4mJ- S )C 9 ° ? », 
^( S ,^)=2 + 4m|^( S ). 

The quark box form factors for = A read 



2s(t + u)C° q ° qq (s) + 2 (t + u + ± j [(t - z)Cf qq (t) + (t- h)Cf qq {t) 



N[t + u + 



A 



+ 2m^A 



2 + 



N 



- 2 



t 



3t - u + —it - uf 



D h qqqq(t, u) - 4 (szh + 7^) t) , 

s(t + u)Cf qq {s) - XC h q z qq {s) 

(t - z)C*(t) + (t - h)CZ(t)] 



2(zN-m 2 q X) D^(t,u) 
3t-u + ^(t-uf 



r ++ 



+ 2 |st 
/ z — h - 



+ 2m 2 n X\D^(s,t), 



Vx 



A; 



qqqq 
1 



(t-^[^ qq (s)--D^(t,u) 



— s 



2 + -(t-u) 



jjhzOO, 



qqqq 



s,t) 



^ l\ h-t 1 
+ 2(t-z)lx z — + 7 = 



2 + 



(t-u)(t-h) 
TV 



Cqqq(t) 



(B.8) 



(B.9) 
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, , f , h — u 1 

_ 2(t _ ft) Az _ + _ 



2 - 



Fi_ = A 



t — u 



-2{t-z){\ z 



s i 

N {S - 
t-h 
N 



1 

7x 



2( f - A )h z — + 7 = 



2 + 
2 - 





N 


,00 / ^ 


- tD hz00 ( 

lLJ qqqq \ 


(t~ 


u){t-h) 




N 




■ u)(u — h) 



N 



(B.10) 
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(b) 

Figure 2: One-loop Feynman diagrams for gg — > Z<p, with (a) <p = h°, H° and (b) <p = A , 
due to virtual quarks and squarks in the MSSM. 
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(a) 

Figure 3: Total cross sections a (in pb) of pp — > Zh° + X via qq annihilation (dashed 
lines) and gg fusion (solid lines) at the LHC (a) as functions of myfi for tan ft = 3 and 30; 
and (b) as functions of tan f3 for m^o = 300 GeV and 600 GeV. For comparison, also the 
Drell-Yan contribution to qq annihilation (dotted lines) is shown. 
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Figure 3 (Continued). 
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Figure 4: Total cross sections a (in pb) of pp — > ZH° + X via qq annihilation (dashed 
lines) and gg fusion (solid lines) at the LHC (a) as functions of m#o for tan/3 = 3 and 
30; and (b) as functions of tan/3 for m^o = 300 GeV and 600 GeV. For comparison, also 
the Drell-Yan contribution to qq annihilation (dotted lines) is shown. 
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Figure 5: Total cross sections a (in pb) of pp — > ZA° + X via bb annihilation (dashed 
lines) and gg fusion (solid lines) at the LHC (a) as functions of m^o for tan (3 = 3 and 30; 
and (b) as functions of tan (3 for m^o = 300 GeV and 600 GeV. For comparison, also the 
quark loop contribution to gg fusion (dotted lines) is shown. 
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Figure 5 (Continued). 
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